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Abstract — Radiometric brightness temperatures below about 12 
GHz provide accurate estimates of path attenuation through 
precipitation and cloud water. Multiple brightness temperature 
measurements at X-band frequencies can be used to estimate rainfall 
rate and parameters of the drop size distribution once correction for 
cloud water attenuation is made. Employing a stratiform storm 
model, calculations of the brightness temperatures at 9.5, 10 and 12 
GHz are used to simulate estimates of path-averaged median mass 
diameter, number concentration and rainfall rate. The results indicate 
that reasonably accurate estimates of rainfall rate and information on 
the drop size distribution can be derived over ocean under low to 
moderate wind speed conditions. 

Index Terms — airborne radiometer, drop size distribution, 
microwave radiometry, rain rate estimation, spaceborne 
radiometer 

I. INTRODUCTION 

For frequencies up to about 12 GHz, the radiometric 
brightness temperature from rain over ocean is closely related 
to the path integrated attenuation (PIA) for low to moderate 
surface wind speeds. Multi-frequency measurements of 
brightness temperature at X-band (8.2-12.4 GHz) over a 
bandwidth of 2-3 GHz, in principle, provide sufficient 
information to estimate parameters of the drop size distribution 
and rainfall rate. The purpose of this paper is to derive 
equations for and investigate the feasibility of this type of 
parameter estimation. 

The instrument concept is similar to that used in the stepped- 
ffequency microwave radiometer (SFMR). In the SFMR, 
multiple bands in the frequency range from 4.6 to 7.2 GHz are 
used to estimate rainfall rate and near-surface wind speed over 
the ocean [1, 2]. Here, the emphasis is different in the sense 
that the objective is to estimate parameters of the path- 
averaged drop size distribution (DSD) and rain rate using X- 
band frequencies. The approach considered here is also 
related to that used for estimation of DSD parameters from 
multi-frequency transmission measurements along a 
microwave link [3-5]. In essence, two path-attenuation 
measurements yield two parameters of the exponential form of 
the DSD or, equivalently, two parameters of the gamma DSD 
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with the shape parameter, p, fixed or expressed as a function 
of one of the variable parameters [6]. There are, however, 
several complicating factors in the application of this 
technique to microwave radiometry. The first is that 
conversion of brightness temperature, T B , to path-integrated 
attenuation. A, is not one-to-one: changes in the characteristics 
of the mixed-phase particles and scattering contributions from 
the ice and water introduce uncertainty in the conversion of Tb 
to A that translate into estimation errors. Cloud liquid water 
presents a somewhat different problem. Cloud droplets are 
Rayleigh scatterers/absorbers within the frequency band of 
interest so that the functional dependence of the specific 
attenuation on frequency is known. This contribution can be 
eliminated in part by considering the difference of path 
attenuations with suitable normalizations. Once the rain 
parameters are estimated from differential quantities, the cloud 
liquid water, in principle, can be recovered from the equation 
for total path attenuation. 

For operation below or above X-band, the error sources 
become larger. Below X-band, the differential attenuation is 
typically small and the relative errors in the estimate render the 
estimates inaccurate. While the differential attenuation is 
relatively strong above 12 GHz, greater scattering 

contributions to T B from ice and snow, as well as the rain 
itself, introduce an increasing amount of variability into the T B 
- A relationship, making the estimation of DSD parameters 
impractical. 

The use of closely-spaced frequencies has also been 
considered for airborne and spaceborne weather radars [7,8]. 
For radar, the differential reflectivity factor serves as an 
estimator of the median mass diameter, D 0 , of rain as well as 
snow. For the X-band radiometer, a ratio of normalized 
differential path attenuations, derived from brightness 
temperatures, provides an estimator for the path-averaged 
median mass diameter of the rain. In a sense, the radiometer- 
based algorithm is more straightforward: for the radar 
application, correction for attenuation must be made before the 
differential reflectivity can be used whereas for the radiometer 
algorithm, a function of the differential attenuation serves as a 
direct estimator of the path-averaged D 0 in rain. As in most 
applications of air- or spaceborne radar and radiometer to 
precipitation, the primary advantage of the radar is its range- 
profiling capability while the attractions of the microwave 
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radiometer are higher reliability, lower cost, and, typically, 
more rapid scanning capabilities. 

II. ALGORITHM CONSIDERATIONS 

The first objective is to express parameters of the DSD and 
rain rate as functions of the PIA or differential PIA. In the 
following section, these quantities are related to the brightness 
temperatures so that a connection is made between the 
measurements and the quantities to be estimated. 

Let A(f) be the PIA (dB) at frequency / (Hz) from the 
storm top (range r=0) to the surface (r = r s ) and let &(/, S ) be 
the specific attenuation (dB km -1 ) at r=s. The quantities 
A(f ) and k{f,s) are related by: 

A(f)=]k(f,s)ds (1) 

0 

The specific attenuation consists of terms corresponding to 
contributions from rain, snow, mixed-phase precipitation as 
well as cloud liquid and cloud ice and various atmospheric 
gases such as water vapor and oxygen. Although the effects of 
atmospheric gases will be assessed in the error analysis given 
later, for the purpose of constructing estimators of the DSD 
parameters, we assume that the contributions from rain and 

mixed-phase precipitation, k pJ and that from cloud liquid 
water, k c , dominate so that 

k = k c +k p (2) 

As the cloud water droplets much smaller than the wavelengths 
of interest, k c can be related to the cloud water content, M c (g 
m- 3 ) by [9] 


where T in (5b) is the estimated mean temperature of the 
cloud droplets. Normalizing (3) by OC n gives: 

K/a n =q n M c ( 6 ) 

where n= 1 or 2 and where 

q l =lm(-K(f,T)) (7a) 

q 2 =[lm(-K(f,T)/kn(-K(f,T))] (7b) 

Dividing A(f ) by (X n and letting A(f)/(X n = A n (f ) then, 

on using (l)-(3), the normalized path attenuation can be 
written 

4 (/) = ][a;\fJ)k p (f,s)^ n (fJ,T)M c (s)Vs ( 8 ) 

0 

If brightness temperature measurements are available at 
frequencies f- > f t ) ) then the difference of the 

normalized path attenuations can be written: 

(9) 

Assuming for either value of n that 

JW„ (/; ,T,T)-q n (f i ,TJW c (s)ds « 0 (10) 

0 

then an approximation for the normalized differential path 
attenuation, independent of cloud liquid water, is: 


k c = [(0.4343 x 67 # )/c] lm(-K)M c (3) 


where c is the speed of light (cm s' 1 ) and K is the dielectric 
factor related to the complex index of refraction of water, m , 
by: 




m 2 -1 

w 2 + 2 


(4) 


The imaginary part of the dielectric factor, Im(— K ) , is a 

function of frequency and temperature, T, or height. To 
remove, in an approximate sense, the influence of cloud water 
on the DSD and rain rate estimates, consider the following two 
normalization factors: 

<*,(/) = [(0.4343 x 6#)/c] (5a) 

a 2 (f,T) = [(0.4343 x6#)/c]Im(— )) (5b) 


<4( (ID 

0 

To understand how information on the raindrop size 
distribution can be obtained from the measurements, we write 
the drop diameter distribution, N(D , s) [m' 3 mm* 1 ], at r=s, as 

N(D,s) = N t (s)n(D,s ) (12) 

where N t is the number concentration (m' 3 ). For the log- 
normal distribution n(D , s) can be expressed as [10]: 

n w {D-,ri, 0 ) = (l/^c£>)exp\ r (lnD-Tj) 2 /2(7 z ] (13) 

and for the Gamma distribution as [1 1]: 
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n c (D;//,A) = [A (fi+1) D M !T(ji + l)]exp[-AZ>] (14) 

where, in general, r|, a for the log-normal distribution and A, p 
for the Gamma distribution are functions of height. In the 
numerical results presented later, we use the median mass 
diameter, D 0 , where [10, 11]: 

A D 0 = 3.67 + JJ. (15) 

for the Gamma distribution and 

D 0 =exp [/7 + 3cr 2 ] (16) 

for the log-normal distribution. 

In the following equations we assume the Gamma 
parameterization; equations for the log-normal 
parameterization can be obtained by replacing n G (D' 9 jU>A) 

with n LN {D\7J i <j ) . The specific attenuation from 
precipitation, k p , can be expressed in terms of the drop 
diameter distribution N(D , s) and the extinction cross section, 
(T e (f,D), by 

k p (f,s) = c k N t (s)jcr e (f,D)n G (D;{i,A,s)dD (17) 

D 

Noting that k p (f,s ) is in dB km" 1 and taking (J e (f,D) to 

be in mm 2 and A, in m" 3 then C k = 4.343x10 3 . Substituting 
(17) into (11) gives 

r s 

&n(fjJi)=C k jNMfc(I*M.^Wn(fjJoV)dlVs 08) 

0 D 

where 

F n (f j ,filD)=a~ 1 (f j ,T)a e (f j ,D)-a n ~ 1 (f i ,T)cr e (f i ;D) (19) 

Although not written as such, F n (fj , f . ; D) varies with height 
because of the temperature dependence of CT e (/, D) . Just as 

it is necessary to evaluate CC n zt some mean cloud droplet 
temperature, a mean raindrop temperature must be used to 
evaluate F n {fj 9 f i \D) in the estimates. Implicit in the 

equations below is the assumption that F n {fj 9 f t \D) is 

independent of height and that, for the estimation procedure, 
the extinction cross sections must be computed at a fixed 
temperature. Errors that arise from this will be discussed in 
section 5. 

An interchange in the order of integration in (18) gives 


0 

K(fjJi)=c k ifj ,/;■;£){ JV, (iX (A /a A s)ds)dD (20) 

D 0 

Replacing the inner integral with the approximation: 

r s 

N t n G (D-,ju, A) = r~ l jN,(s)n G (D;fi, A,s)ds (21) 

0 

allows (18) to be written: 

Anifj.fi) = c k r s N t jtt G (D;/i,~K)F n (f j ,f i ;D)dD (22) 

Although N t and A yield the path-averaged number 

concentration and slope parameter only for range-independent 
drop size distributions, these will be taken to be the estimated 
values of these quantities. 

If we have access to brightness temperatures measured at 3 
frequencies, then, using the approximation (22), the ratio 

)(<*«. A)) 

yields a quantity independent of the mean number 
concentration and approximately independent of cloud water: 

AnUjjA} _ ljD)dD ^ 

A n if, 2 . /, ) " jn G (D;ji,~K)F n (/ 2 , f x ; D)dD 

If (X is a constant or a function of A , then (23) can be solved 
numerically for A or D 0 . Once A has been found, N t follows 

from (22). Next, k p can be estimated from N t and A : this 
determines the first term of (8). Estimating A(f ) from 
T B (f ) and using estimates of N t and A from (22) and (23), 

respectively, provides the integrated cloud water content, I M . 
Specifically, from (8), with q n in (7) taken to be unity, an 
estimator of I M is: 

I M = \m c (s)ds =a n ~ l (f,T r )[A(f) — 

o (24) 

c k r s N, §n G (D,ju,A)(T e (D, f)dD] 

D 


where f can be any of the frequencies at which the brightness 
temperature is measured. 
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An estimate of the path-averaged rain rate (mm h" 1 ) follows 
from the definition of this quantity along 

with N t and A obtained from the previous equations: 

R = c R N t jn G (D;u,A)v(D)D 3 dD (25) 

D 

where the velocity distribution of raindrops, v(D ) (m s' 1 ) is 
approximated by [12]: 

v(D) = 9.25[1 - exp(-0.068D 2 - 0.488Z))] (26) 

where D is in mm in (25) and (26) so that 

C R = 0 . 62 TXlCr 3 (27) 

As shown by the equations above, measurements of path 
attenuation are needed at a minimum of 3 frequencies to 
estimate the path-averaged rain rate and the parameters 

N t and A of the drop size distribution. Although (23), which 
provides an estimate of D Qj is independent of the effective 
range, r s , through the precipitation, this is not the case for 

N t and R . In the case of stratiform rain, the effective range 

would include the melting layer and rain but not the dry snow 
above the melting layer, an estimate of which would require 
either the detection of the melting layer by radar or an estimate 
of the surface temperature and lapse rate. 

III. Conversion of Brightness Temperature to Path Attenuation 

To use the equations given in the previous section, the 
brightness temperatures must be related to the PIA or 
differential PIA. We begin by relating the brightness 
temperature, T B (/) , to the PIA by [13, 14]: 

Mf) = c 0 (f) + Cl (f)ln[T 0 -T B (f)] . (28) 

where T () is a constant brightness temperature. Recalling that 
the differential normalized PIA, ) , is defined by 

(9), and letting 

y = &K(fjJi) (29) 

then (28) implies that y can be expressed in the following form 

y = £ 0 +£ i*i+&*2 (30) 


*i=ln lT 0 -T B (fj)] (31) 

x 2 =ln[T Q -T B (fi)] (32) 

Although the £ coefficients in (30) can be expressed in terms 
of c 0 and ci in (28), a more accurate procedure is to generate 
A(f i ),T B (f i );i ' = 1,2,3 from model storms and determine 
the £ and c coefficients by linear regressions. 

IV. Storm Model and Brightness Temperature Calculations 

We construct a simple model of stratiform rain by the 
following procedure [6]. From ground-based disdrometer- 
measured drop size distributions the number concentration and 
median mass drop diameter, D 0 , are calculated. Next, a best- 
fit p is computed either for each DSD or for the ensemble. For 
the results here, we use p=2. As can be seen from (14) and 
(15), this information specifies a gamma distribution of drop 
diameters which is then used to characterize the distribution 
along an entire vertical column of precipitation. The 
procedure is repeated for a series of measured drop size 
distributions so that a set of such gamma distributions (and the 
corresponding vertical columns) are generated. The height of 
the rain layer is fixed at 4 km; above this, a melting layer and 
snow layer are appended. In the melting region, the model of 
Yokoyama and Tanaka [15] is used along with the effective 
medium or Bruggeman approximation [16] for the effective 
dielectric constant of mixed-phase hydrometeors. Above the 
0° isotherm a 1 km layer of snow is added where the particle 
mass density is taken to be either 0.05 g cm' 3 or 0.2 g cm' 3 . 
Throughout the vertical column the mass flux, or equivalent 
rain rate, is taken to be constant. Added to the precipitation- 
sized particles is cloud liquid water which is specified by the 
liquid water content (g m" 3 ) as a function of height. For the 
results presented here, a uniform distribution of cloud water of 
1 km depth is located above the 0° isotherm. M c is allowed to 

take on values from 0 to 2 g m' 3 so that the integrated cloud 
water can range from 0 to 2 kg m' 2 . The relative humidity is 
allowed to take on values between 70% and 100%. In the 
radiative transfer model of Kummerow [17], the surface 
emissivity is computed as a function of wind speed and 
frequency using the models of Hollinger et al. [18] and 
Stogym [19]. In the baseline model, we take the wind speed to 
be 7 m s' 1 ; the effects of changes from this value are discussed 
later. 

From the storm and surface scattering models, brightness 
temperatures and total path attenuations at various frequencies 
are calculated over each vertical column. The brightness 
temperatures are calculated using the Eddington 
approximation for the radiative transfer equation [17]. 
Although arbitrary incidence angles can be used, we restrict 
the calculations to nadir incidence. 


where 
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V. Results 

Using the storm model just described, we calculate the 
apparent radar reflectivity factor profiles, 
dBZ m = 101og IO Z m , at 9.5 GHz for nadir incidence above 
the storm where the apparent, Z m , and actual radar reflectivity 
factor, Z , are related at radar range r by 

r 

Z m (r) = Z(r)exp[-0.21nl0 fk(s)ds] • (33) 

0 

A sequence of 379 profiles is shown in the top panel of Fig. 1. 
In this case, the integrated cloud water is taken to be 0.5 kg m' 2 
with snow density of 0.2 g cm' 3 and 80% relative humidity. 
The radar bright-band, corresponding to mixed-phase 
hydrometeors in the melting layer, is evident in the region just 
below 4 km. In the second panel from the top, the reflectivity 
factor difference, 

(12,9.5 GHz) = dBZ m (l2GHz) - dBZ m (9.5GHz) (34) 

is shown. The characteristics of the reflectivity and 
differential reflectivity profiles are discussed in [6]. The 
corresponding brightness temperature results are shown in the 
lower two panels: T B (/ = 9.5GHz) 

and ST b (1 2GHz,9.5GHz) = T b (Y2jGHz)-T b {95GHz ) . 

Notice that 8T b is bounded below at approximately 6 K. The 

bulk of this difference is caused by the differential path 
attenuation from cloud liquid water, water vapor and 
molecular oxygen. If the integrated cloud water is set to zero 
the minimum brightness temperature difference is about 4 K; if 
it is set to 1 kg m' 2 , the minimum 8T b is about 8 K. In other 
words, at light rain rates, an increase in the integrated cloud 
water of 1 kg m' 2 results in an increase in the differential 
brightness temperature between 9.5 and 12 GHz of about 4 K. 

To estimate N t , D 0 and R from (22), (23), and (25) it is first 

necessary to convert the T B to A and using the functional 
forms given by (28) and (30) for the storm model data shown 
in Fig. 1. For all fits, we assume that T 0 = 280 K. Tables 1 

and 2 provide the coefficients Cq, Ci and £i, and 2^, 
respectively for several frequencies and frequency pairs. Plots 
of path attenuation versus brightness temperature are shown in 
Fig. 2 for the frequencies 9.5, 10 and 12 GHz along with the 
fits (solid lines) using the coefficients from Table 1. It should 
be pointed out that the solution set of (30) defines a plane in y- 
x r x 2 space: the solid curve shown in Fig. 3 is obtained by first 
expressing T B (f = 12G/fe)as a linear function of 
T B (f = 9.5 GHz) and then plotting y in (30) as a function of 
T B (f = 9.5 GHz) . It should also be noted that in Fig. 3 and 
for the coefficients in Table 2 we have chosen the first 
normalization factor so that S\ n = . This is also the 

normalization used for the D 0 estimates. 



250 r 


200 


150 


100 






Fig. 1 : Vertical profiles of measured radar reflectivity factor, 
Z m (9.5 GHz) (top); profiles of differential radar reflectivity 
factor, Z m (12 GHz, 9.5 GHz), (second from top); brightness 
temperature, T B ( 9.5 GHz) (third from top), and differential 
brightness temperature, T B (12 GHz)- T B (9.5 GHz) (bottom) 
for the baseline storm and surface model. 



Fig. 2: Path-integrated attenuation (unnormalized) (dB) at 3 
frequencies versus corresponding brightness temperature, T B , 
using the baseline stratiform storm model. 
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Fig. 3: &\ x (12, 9.5 GHz) versus T B (12 GHz) and T B (9.5 

GHz) where values are derived from model data shown in Fig. 
1 ; solid line represents a fitted curve through the data. 


Table 1: Coefficients c 0 (f) , c, (/) in 

Mf) = C 0 (f ) + c, C f)ln[T 0 - T b (/)] ; T 0 =280 K 


Frequency (GHz) 

Co 

Cl 

9.5 

11.89 

-2.33 

10 

11.99 

-2.35 

11 

12.25 

-2.40 

12 

12.63 

-2.48 


Table 2: £ coefficients, with T 0 

=280 K 



4o + MT 0 

r„(/,)+£in(r 0 

(fpfi) GHz 


#1 

£2 

(10, 9.5) 

0.019 

-0.137 

0.133 

(11,9.5) 

0.093 

-0.115 

0.096 

(12, 9.5) 

0.216 

-0.094 

0.051 

(11,10) 

0.065 

-0.104 

0.091 

(12, 10) 

0.168 

-0.085 

0.052 

(12,11) 

0.085 

-0.063 

0.046 




Do 

Fig. 4: Ratio of frequency-normalized differential path- 

attenuations (top) and corresponding estimate of this quantity 
derived from brightness temperatures (bottom). Results are 
plotted as a function of the median mass diameter, D 0 . The 
solid line represents the best fit curve of D 0 as a function of 
the ratio plotted along the ordinate. 

The coefficients given in Tables 1 and 2 are fixed for all 
calculations. Although these are the best-fit coefficients for 
the baseline model results, they are not the best-fit coefficients 
for models with changes in cloud water, water vapor, snow 
density or surface wind speed. Clearly, the feasibility of the 
method depends not only on the quality of estimates using data 
from a particular model but the accuracy of estimates from 
storm models that differ from the baseline. 

Using (23), we compute and plot in Fig. 4 the ratio 
p = a,(/3./,)/<^(/ 2 ./,) versus D 0 where 
(/i,/ 2 ,/ 3 ) = (9.5,10,12 )GHz . For the scatter plot in the 

top panel, the true differential normalized path integrated 
attenuation is used to compute the ratio while for the bottom 
panel S\(f 3 , f r ) and &\(f 39 f 2 ) are estimated from the 

brightness temperatures using (30) with the coefficients listed 
in Table 2. Differences between the upper and lower scatter 
plots are largely due to errors incurred in converting brightness 

temperatures to Notice also that at small D 0 , (typically 
light rain rates), the true value of p is determined primarily 
by the differential attenuations of the cloud and atmospheric 
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gases and is nearly independent of the characteristics of the 
precipitation. Taking the data in the bottom panel of Fig. 4 as 

the basis of estimating D 0 from p , then on using a quadratic 
fit to the data, we obtain: 

D 0 =17.02-4.81p + 0.350/> 2 (35) 

where (35) is plotted as a solid line in the lower panel of Fig. 
4. 

The results shown in Figs. 1-4 address the forward problem 
where the brightness temperature and path-integrated 
attenuation is computed for a given frequency and storm 
model. The remainder of the paper is focused on the inverse 
problem where estimates of rain rate and parameters of the 
drop size distribution are obtained from brightness 
temperatures at several frequencies. 

To obtain an estimate of Z) 0 , we first convert the T B to 
(/ 3 , fi ) and (/ 3 , f 2 ) using (30) and the coefficients in 

Table 2. From these quantities p is computed. D 0 then 
follows from (35). The results of the procedure are shown in 
Fig. 5 (top). Once D 0 is obtained, N t is determined from 
(22). A scatter plot of the logarithm (base 10) of the estimated 
N t values versus the ‘true’ values (i.e., those values used in 
the forward calculations of T B and A) is shown in the bottom 
panel of Fig. 5. Less scatter in the N t estimates is achieved by 

using SA 2 rather than S\ so that the n=2 normalization in 
(22) has been used. It worth emphasizing that although (35) is 
the best-fit curve of A, versus p for the data generated from 

the baseline model only, the coefficients of the fit are used for 
all subsequent estimates despite changes in the storm model 

parameters and wind speed used in the calculations of T B . 

From the estimated values of Z) Q and N t , shown in Fig. 5, the 

rain rate estimates are calculated from (25). Plots of the 
results versus sequence number are shown in the top panel of 
Fig.6; below this are plotted the true rain rates. A scatter plot 
of the estimated versus true rain rates is given in the bottom 
panel of Fig. 6. 

As mentioned in section 3, evaluation of the integrals in (22) - 
(25) requires computations of CT,(/,D) and F n {f 

at a fixed temperature. For the results in Figs. 4 and 5, the 
near-surface temperature of 24 C was used. If we change this 
to the temperature of the midpoint of the rain layer (12 C), we 

obtain almost identical results for D 0 but significantly 

different values for N t &ndR . Plots of the estimated versus 
true values of these latter two variables are shown in Fig. 7; it 



Fig. 5: Top: Estimated median mass diameter (D 0 ) versus the 
true value using the fitting curve shown in Fig. 4. Bottom: 
logarithm of estimated number concentration, N t ,. 



O 95 190 284 379 



Rain Rate (mm hr ’) 


Fig. 6: Top: Estimated and true rain rates versus sequence 
number. Bottom: Scatter plot of estimated and true rain rates 
from data in the top panels. 
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can be seen that a change in the assumed temperature of 
raindrops from 24 to 12 C increases the estimates of number 
concentration and produces positively biased rain rates. For 

subsequent plots and computations, <J e (f,D) and 

F n ( f j , f t ; D ) are evaluated at 24 C. 



0 5 10 15 

Rain Rate (mm hr -1 ) 


Fig. 7: Effects of a change in temperature from 24 C to 12 C at 
which estimates are evaluated. Top: Estimated versus true 
values of log(Nt); Bottom: estimated versus true values of rain 
rate. 

Before discussing changes to the baseline model, we note that 
estimates of integrated cloud water using (24) are generally 
unreliable. The reason for this can be understood by 

recognizing that (24) expresses the cloud water 

attenuation, 0C c (f 9 T c )I M , as a difference between the total 

attenuation and the attenuation from precipitation alone. The 
estimate is positively biased because attenuation contributions 
from water vapor and 0 2 have not been subtracted from the 
total. This problem can be partly circumvented by modifying 
the estimate to include this subtraction using nominal values 
for these contributions. However,, a more serious problem 
with (24) is that for moderate and heavy rain rates, the cloud 
water attenuation, for the models considered here, is a small 
fraction of the total. As a consequence, variations in the 
estimate of the second term on the right-hand side of (24) 


introduce large errors in I M because of the subtraction of 
quantities of roughly equal magnitude. 



Fig. 8: Estimated versus true values of D 0 , log(N t ) and R for a 
modified stratiform storm model: snow density= 0.2 g cm' 3 , 
integrated cloud water = 0.25 kg m 2 , RH=90%. 



Fig. 9: Estimated versus true values of D 0 , log(Nt), and R for a 
modified stratiform storm model: snow density= 0.05 g cm' 3 , 
integrated cloud water = 0.75 kg m' 2 , RH=90%. 
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The remaining examples focus on results from modifications 
of the baseline stratiform storm model. Shown in Figs. 8-10 

are scatter plots of the estimated versus true values of Z) Q 

(top), log( N t ) (center) and R (bottom) for the sets of model 

parameters in Table 3. (It is worth noting that a replacement 
of the effective medium with the Maxwell Garnet formulation, 
water matrix with ice inclusions, for melting snow particles 
produces a significant increase in brightness temperatures 
which leads to increases in N t and positive biases in the rain 
rate. The results are somewhat similar to those shown in Fig. 
7. A detailed study of the effects of the mixed-phase region is 
beyond the scope of the paper, however.) 

Table 3: Parameters used in the calculation of brightness 
temperature 


Model 

Snow 
density 
(g cm' 3 ) 

Integrated 
cloud 
water 
(kg m' 2 ) 

Relative 

humidity 

(%) 

Ocean 
wind 
speed 
(m s' 1 ) 

Relevant 

figures 

baseline 

0.2 

0.5 

80 

7 

1-7 

Modi 

0.2 

0.25 

90 

7 

8 

Mod 2 

0.05 

0.75 

90 

7 

9 

Mod 3 

0.2 

2 

80 

7 

10 

Mod 4 

0.2 

0.5 

80 

15 

11 



Fig. 10: Estimated versus true values of D 0 , log(N t ), and R for 
a modified stratiform storm model: snow density= 0.2 g cm' 3 , 
integrated cloud water = 2 kg m~ 2 , RH=80%. 


Despite the fairly large variability in D 0 and N t in these 

examples, the rain rates tend to have a small amount of scatter 
relative to the true values with generally negative biases for 
cases shown in Figs. 8 and 9 and a positive bias for the case in 
Fig. 10. As noted earlier, the bias is also affected by the 
assumed temperature at which CT e (/,£>) and 

F n {f j , fi ; D) are evaluated and the type of normalization 

used to form the estimates. The reason for the relatively 
accurate rain rate estimates appears to arise from the fact that 

D 0 and N t are inversely correlated so that an overestimate in 

one is compensated by an underestimate in the other. This 
compensation works well for the rainfall estimate because of 
the underlying path attenuation constraint; that is, the 
brightness temperatures at the various frequencies are highly 
correlated with attenuation. Moreover, because rain rate and 
attenuation are approximately equal to the same moment of the 
drop size distribution [20], this constraint extends, in an 



O 95 190 284 379 



Fig 11: Estimated versus true rain rates for an ocean wind 
speed of 15 m s' 1 . 

approximate sense, to the rain rate estimates. For similar 
reasons, a direct conversion of a brightness temperature 
measurement at X-band to path-averaged rainfall rate is fairly 
accurate. The primary advantage of the use of multiple 
frequencies is the added information provided on path- 
averaged number concentration and median mass diameter. 
Although the rainfall estimates appear to be reasonably 
accurate, the improvement over a straightforward rain rate 
estimator using a single brightness temperature needs to be 
explored. 

In previous examples, the T B were calculated under the 

assumption of an ocean wind speed of 7 m s' 1 . In the baseline 
model, the assumed wind speed was also taken to be 7 m s' 1 . If 








TGRS-001 5 1-2004 


10 


the T b are calculated at a higher wind speed than the assumed 
value of 7 m s' 1 , we find that D 0 is positively biased and N t is 

negatively biased. Fig. 11 shows the effect on the rain rate 
estimates for a wind speed of 15 m s 1 , where, as always, a 7 m 
s' 1 speed is assumed in the retrieval. Note that the relative 
error is particularly large at light rain rates. For low winds, the 
errors in the estimated quantities are much smaller. In general, 
the results suggest that errors are acceptable up to wind speeds 
of about 10 m s' 1 to 15 m s' 1 although accuracy depends on the 
magnitude of the rain rate and cloud liquid water as well as 
other error sources. Some improvement in the retrieval should 
be possible by estimating wind speed in nearby rain-free areas. 
However, in high-wind environments such as hurricanes, the 
algorithms are unstable and the technique devised for the 
SFMR becomes applicable [1,2]. 

VI. Summary and Conclusions 

Brightness temperature measurements at multiple frequencies 
within X-band, over a span of 2-3 GHz, may provide 
information on the path-averaged rainfall rate and parameters 
of the raindrop size distribution. The potential for this type of 
estimation arises from the fact that at X-band, brightness 
temperature is well correlated with path-integrated attenuation 
and that a difference of path attenuations normalized by 
frequency is nearly independent of cloud liquid water 
attenuation. Moreover, a ratio of such differences, derived 
from brightness temperatures at three X-band frequencies, 
provides an estimate of the median mass diameter, D 0 , of the 
rain. The equations also lead to estimates of the mean number 
concentration, N t , and rainfall rate. A simple simulation using 
a stratiform storm model suggests that reasonably accurate rain 
rates are possible even though the variability in D 0 and N t can 
be high. On the other hand, estimates of integrated cloud 
water content appear to be unreliable. Although the approach 
appears to be feasible for low to moderate wind speeds over 
ocean, more detailed studies are needed on variations in the 
shape parameter, jlx, changes in the drop size distribution with 
height and the effects of the melting layer. The error budget 
must also include the measurement error of the brightness 
temperatures in the context of the system design, measurement 
requirements and scanning strategy. 

ACKNOWLEDGMENT 

I wish to thank Chris Kummerow for the use of his radiative 
transfer code and Jim Weinman for helpful discussions and 
suggestions. 

REFERENCES 

[1] P.G. Black and CL. Swift, “Airborne stepped frequency microwave 
radiometer measurements of rainfall rate and surface wind speed in 
hurricanes,” 22 nd Conf. Radar Meteorology, pp. 433-438, Sept., 1984, Amer. 
Meteor. Soc., Boston, MA. 

[2] E.W. Uhlhom and P.G. Black, ‘Verification of remotely sensed sea 
surface winds in hurricanes,” J. Atmos. Oceanic TechnoL, vol. 20, pp. 99- 
116, Jan. 2003. 


[3] Y.T. Furuhama and T. Ihara, “Remote sensing of path-average raindrop 
size distributions from microwave scattering measurements,” IEEE Trans. 
Antennas Propagate vol. AP-29, pp. 275-281, Mar. 1981. 

[4] T. Ihara, Y. Furuhama, and T. Manabe, “Inference of raindrop size 
distribution from rain attenuation statistics at 12, 35, and 82 GHz,” Trans. 
IECE Japan, vol. E67, pp. 211-217, April 1984. 

[5] R.F. Rincon and R.H. Lang, “Microwave link dual-wavelength 
measurements of path-average attenuation for the estimation of drop size 
distribution and rainfall,” IEEE Trans. Geosci . Remote Sens., vol. GE-40, pp. 
760-770, April 2002. 

[6] G. Zhang, J. Vivekanandan, and E. Brandes, “A method for estimating 
rain rate and drop size distribution from polarimetric radar measurements,” 
IEEE Trans . Geosci. Remote Sens., vol. 39, pp. 830-841. 

[7] R. Meneghini, L. Liao, S. W. Bidwell, and G. M. Heymsfield, “On the 
feasibility of a Doppler weather radar for estimates of drop size distribution 
using two closely-spaced frequencies,” IEEE Trans. Geosci. Remote Sens., 
vol. GE-39, Oct. 2001. 

[8] R. Meneghini, S.W. Bidwell, L. Liao, R. Rincon, and G.M. Heymsfield, 
“Differential-frequency Doppler weather radar: theory and experiment,” 
Radio Science , vol. 38, pp. 5.1-5.10, Mar. 2003. 

[9] L.J. Battan, Radar Observations of the Atmosphere. The Univ. of 
Chicago Press, 1973, 324 pp. 

[10] G. Feingold and Z. Levin, “The log-normal fit to raindrop spectra from 
frontal convective clouds in Israel,” J. Climate Appl. Meteor., vol. 25, pp. 
1346-1363, 1986. 

[11] C.W. Ulbrich, “Natural variations in the analytical form of the raindrop 
size distribution,” J. Climate Appl. Meteor., vol. 22, pp. 1764-1775, 1983. 

[12] R. Lhermitte, “Cloud and precipitation sensing at 94 GHz,” IEEE 
Geosci. Remote Sens., vol. GE-26, pp. 207-216, Jan. 1988. 

[13] J.A. Weinman, R. Meneghini, and K. Nakamura, “Retrieval of 
precipitation profiles from airborne radar and passive radiometer 
measurements: comparison with dual-frequency radar measurements,” J. 
Appl. Meteor vol. 29, pp. 981-993, Oct. 1990. 

[14] Durden et al., “Measurement of rainfall path attenuation near nadir: a 
comparison of radar and radiometer methods at 13.8 GHz,” Radio Science , 
vol. 30, pp. 943-947, July- Aug. 1994. 

[15] T. Yokoyama and H. Tanaka, “Microphysical processes of melting 
snowflakes detected by two-wavelength radar,” J. Meteor. Soc. Japan , vol. 
62, pp. 650-666, 1984. 

[16] D.A.G. Bruggeman, “Berechnung vershiedener physikalischer 
Konstanten von heteorogenen Substanzen. I. Dielectrizitatskonstanten und 
Leitfahigkeiten der Mischkorper aus isotropen Substanzen,” Ann. Phys., vol. 
24, pp. 636-679, 1935. 

[17] C. Kummerow, “On the accuracy of the Eddington Approximation for 
radiative transfer in the microwave frequencies,” J. Geophys. Res., vol. 98, 
pp. 2757-2765, Feb. 1993. 

[18] J. Hollinger, R. Lo, G. Poe, R. Savage, and J. Pierce, Special Sensor 
Microwave Imager User’s Guide , 120pp., Naval Research Laboratory, 
Washington, D.C., 1987. 

[19] A. Stogym, “The emissivity of sea foam at microwave frequencies,” J. 
Geophys. Res., vol. 77, pp. 1658-1666, 1972. 

[20] D. Atlas and C.W. Ulbrich, “Path- and area-integrated rainfall 
measurement by microwave attenuation in the 1-3 cm band,” J. Appl. 
Meteor., vol. 16, 1322-1331, 1977. 



